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Abstract. We present a method of lensing reconstruction on the full sky, by extending the 
optimal quadratic estimator proposed by Okamoto & Hu (2003) to the case including the 
curl mode of deflection angle. The curl mode is induced by the vector and tensor metric 
perturbations, and the reconstruction of the curl mode would be a powerful tool to not 
only check systematics in the estimated gradient mode but also probe any vector and tensor 
sources. We find that the gradient and curl modes can be reconstructed separately, thanks 
to the distinctive feature in the parity symmetry between the gradient and curl modes. We 
compare our estimator with the flat-sky estimator proposed by Cooray et al (2005). Based 
on the new formalism, the expected signal-to-noise ratio of the curl mode produced by the 
primordial gravitational-waves and a specific model of cosmic strings are estimated, and 
prospects for future observations are discussed. 
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1 Introduction 

Ongoing, upcoming and next-generation experiments of CMB, such as PLANCK [1], PO- 
LARBEAR [2], ACTPol [3], SPTPol [4], CMBPol [5], and COrE [6] would resolve not only 
the fine structure of the primary CMB anisotropies but also tiny effects that are important at 
small scales. One of the most important signals in those effects is the weak lensing: the deflec- 
tion of CMB photons coming from the last-scattering surface by metric perturbations along 
our line-of-sight. Recent studies show that the lensing fields involved in the CMB anisotropies 
are reconstructed from the cross-correlations between CMB and large-scale structure [7, 8], 
and CMB maps alone [9-11], and the lensing information from upcoming experiments would 
provide us an opportunity to probe the late-time evolution of the structure in the Universe 
(e.g., [12-19]). 

Several studies have investigated a method to reconstruct the deflection angle, which 
characterizes the effect of weak lensing on CMB maps (e.g., [20-28]). In general, the deflection 
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angle in a direction n, d(n), where n is the unit vector defined on the unit sphere, is 
decomposed into gradient and curl part as (e.g., [26]) 

d{h) = V(l){n) + (*V)ro(n) , (1.1) 

where the first term, W(f){ii), and second term, (*V)tz7(n), represent gradient and curl mode 
of deflection angle, respectively, and, in the polar coordinate, the covariant derivative on 
the unit sphere is given by V = eg{d/d9) + {sy^/ sm6){d/dip) with eg and describing 
the basis vectors in the polar coordinate. The symbol, denotes a operation which rotates 
the angle of two-dimensional vector counterclockwise by 90-degree; for a vector on the unit 
sphere expressed in terms of the basis vectors, a = ageg + a^e^p, the operator, act on 
a as (*a) = age^p — a^peg [26]. Hereafter, we call the potentials, (p and ro, "scalar lensing 
potential" and "pseudo-scalar lensing potential" , respectively. 

The scalar metric perturbations such as the matter density fluctuations at linear order 
produce only the gradient mode, and the curl mode is usually neglected in the algorithm of 
lensing reconstruction. However, the curl mode can be induced by vector and/or tensor metric 
perturbations. In the conformal Newton gauge, the line element in the polar coordinate 
system is given as 

ds^ = a'^{r]){-{l + 2A)drf - 2Bidr]dx' + [(1 + 2C)7ij + 2Dij]dx'dx^} , (1.2) 

where a is the scale factor, rj is the conformal time, A and C are the scalar components, Bi is 
the vector component (i?j|j = 0), and Dij is the tensor component (-Djj|j = 0, and Da = 0). 
The unperturbed spatial metric, ^ij, is given by 

-fijdx'dx^ = dx^ + X^{dO^ + sin^ 6dip^) . (1.3) 

Then, the pseudo-scalar lensing potential is described by [29] 

tn(n) = (^V)-2 r 

Jo X^smO 

where the quantities, x Xs, are the comoving distance and comoving distance at the last 
scattering surface, respectively, tjq is the conformal time today, and the quantities, Qg and 
^Ip, are defined as 

^aiv,Xn) = Ba{r],xn) + 2D^aiv,Xn) {a = 9,Lp). (1.5) 

According to the above equations, the primordial gravitational- waves produce the curl mode 
(e.g., [30]). The cosmic strings can also produce the curl mode through the vector and tensor 
perturbations (e.g. ,[29, 31]). This implies that the curl mode is a smoking gun of cosmic 
strings and other vector or tensor sources. Even if these sources are absent, the higher-order 
density perturbations and foreground contaminations generate not only the gradient mode 
but also the curl mode. These sources would cause systematics in the estimation of the 
gradient mode. Thus, the evaluation of the contaminations of these sources in the curl mode 
would be helpful to estimate the contributions of these sources in the gradient mode. 

In this paper, we present a method to reconstruct both the gradient and curl modes 
of deflection angle. The reconstruction of the curl mode has been previously discussed in 
Refs.[26, 32]. Ref.[26] proposed an algorithm based on the likelihood analysis. Including 



d^eiVo - X, Xn) dVtp{rjo - x, xn) 
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polarizations, this estimator, in principle, suppresses the noise contribution for the recon- 
structed potentials, but numerically cost, compared to the estimator in Ref.[24]. On the 
other hand, in Ref.[32], they consider the flat-sky limit and empirically define a quadratic 
estimator of the curl mode. Since current and future CMB missions will cover nearly entire 
sky, a full-sky algorithm including the curl mode are highly desirable. In this paper, we derive 
a full-sky estimator of the gradient and curl modes on the full sky, extending the full-sky 
formalism for the gradient mode in Ref.[25]. Then, we compare the fiat-sky estimator with 
the one on the full sky, and show that the empirically defined estimator in Ref.[32] can be 
derived from the full-sky estimator. In addition, based on our full-sky estimator, possible 
implications to detection of curl mode from primordial gravitational-waves and cosmic strings 
are discussed. 

This paper is organized as follows. In section 2, we briefly summarize the lensing effect 
on the CMB anisotropics. In section 3, we extend the quadratic estimator of Ref.[25] to the 
case including both the gradient and curl modes of deflection angle. In section 4, we compute 
the noise spectrum of the estimator in both the full- and flat-sky cases. In section 5, we discuss 
implications for primordial gravitational-waves and cosmic strings by reconstructing the curl 
mode. Section 6 is devoted to summary and conclusion. 

In this paper we adopt the cosmological parameters assuming a flat Lambda-CDM 
model consistent with the results obtained from Ref.[33]; the density parameter of baryon 
n^h'^ = 0.022, of matter n,^h^ = 0.13, dark energy density Oa = 0.72, scalar spectral index 
rig = 0.96, scalar amplitude Ag = 2.4 x 10~^ and the optical depth, r = 0.086. In Table. 1, 
we summarize the meaning and deflnition of the quantities used to reconstruct the lensing 
potentials from observed CMB maps. 



2 Weak lensing of the CMB 

Here we briefly review the lensing effect on the CMB anisotropics in the full-sky case, includ- 
ing both the gradient and curl modes of deflection angle. We flrst discuss the lensing effect 
on CMB temperature in section 2.1, and the similar discussion for polarizations is given in 
section 2.2. The detailed calculation of lensing effect is presented in, e.g., Ref.[34] in the 
absence of the curl mode, and in Ref. [30] including the curl mode. 

2.1 Lensing effect on CMB anisotropies: temperature 

Let us first discuss the lensing effect on the temperature anisotropies. The lensed temperature 
fluctuations in a direction n, Q{n), are transformed into the harmonic space according to 

e,,™ = J dnoyl^{n)G{n) , (2.1) 

with the quantities, Qi^m and o3^f,m.(^)) describing the harmonic coefficients, and the spin- 
spherical harmonics, respectively. The lensed temperature fiuctuations are related to the 
unlensed temperature fluctuations, 0(ri), through 0(n) = Q{h+d), where d is the deflection 
angle. Usually, the deflection angle is a small perturbed quantity, \d\ <^ 1, and the lensed 
temperature fluctuations may be expressed as 

G(n) = e(n) + d- VG(n) + 0(|d|2), (2.2) 
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Table 1. Notations for quantities used to reconstruct the lensing potentials from observed CMB maps. 
The symbols are divided into three categories. The quantities in the middle eight rows are needed 
to compute the optimal quadratic estimator. In the bottom three rows, we describe the quantities 
needed to compute the optimal combination. 



Symbol 


Definition 


Meaning 


Full sky / Flat sky 


Full sky / Flat sky 




X (= (j), w) 


Eq.(l.l) 


Scalar or pseudo-scalar lensing potential 


a (or /3) 




A pair of two CMB maps, X and Y 


/ \ 7 ^ 

/ ^(") 
i,L,L' 1 ^ f.,L,L' 


Eq.(3.10) / (3.42) 


Estimator 


Eq.(3.27) / (3.43) 


Weight function 




Eq.(3.26) / (3.41) 


Noise spectrum 


x,(ol) I 


Eq.(3.21) / (3.40) 




rx,(a) 1 -eX,{a) 
h,L,L' 1 J t,L,L' 


Table.2 / 3 






Eq.(2.7) 




QX 


Eq.(2.13) 




and qS^ 


Eq.(2.16) 




^(c) / ^(c) 


Eq.(3.29) / (3.45) 


Optimal combination 


^x,(c) 1 ^x,(c) 


Eq.(3.30) / (3.46) 


Noise spectrum for optimal combination 


j^x,{a,l3) 1 j^x,(a,l3) 


Eq.(3.31) / (3.47) 


Noise cross-spectrum 



where we expand 0(n -|- d) in terms of the deflection angle d. Hereafter we neglect the 
contributions of 0(|dp) in Eq.(2.2). The harmonic coefficients of the lensed quantities are 
obtained by transforming Eq.(2.2) into the harmonic space, according to Eq.(2.1). Using 
the expression of deflection angle (1.1), the lensed temperature anisotropics in the harmonic 
space are given by [30] 

@L,M = @L,M + J dnoK,A/(^)[V0(n) + {-kV)w{n)] • Ve(n) 

= e^,M + EE0^>'(-l)''(_M,lii') E oSl,,,xe,m, (2.3) 

Am l',m' x=(t>,T^ 

where the quantities, O'S'^^^i and qS^^^', are defined as 

(-M L i) = / ^^oK,M(n)[Vo3^.,™(n)] . [Vo3^,',^'(n)] , (2.4) 

(-1)'' (_M m m') = / ^^oH,A/(^)[(*V)o3^f,.„(n)] • [V o3^,'v(^)] • (2.5) 
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As shown in appendix A (Eqs.(A.7) and (A. 8)), the quantities, O'Sf^^^, and qS^^^,, are 
expressed in terms of the Wigner-3j symbols, and the results are 



o^LAi' = V ["^(^ + 1) + ^(^ + 1) + ^ + 1)] ( ) > (2-6) 



'(2£ + l)(2^' + l)(2L + l 

0-11/ 1 1 -1 



We note that the quantities ^ £/ and oS^^ satisfy 



(2.7) 



The above equations come from the parity symmetry of 0, (j) and ro; the temperature 
anisotropies and the scalar lensing potential are even parity, while the pseudo-scalar lensing 
potential is odd parity. In fact, Eq.(2.8) is checked by changing the variable, n — t- — n in 
the right-hand side of Eqs.(2.4) and (2.5). Under this transformation, the spin-0 spheri- 
cal harmonics are multiplied by a factor (—1)^, and the derivatives become V — t- V and 
(*V) — )■ — (*V), respectively. As a result, the right-hand sides of Eq.(2.4) and Eq.(2.5) are 
multiplied by a factor of (—1)^+^+^ and —(—1)^+^+^ , respectively. Eq.(2.8) is also checked 
with the formulas of the Wigner-3j symbols (see Eq.(A.2) in appendix A). 

From Eq.(2.8), 0*5^ becomes zero if L -|- £ + / is an odd integer, and the coefficient 
qS^^ fi vanishes when is an even integer. These properties are essential for a separate 

reconstruction of gradient and curl modes in subsequent analysis. 

2.2 Lensing effect on CMB anisotropies: polarizations 

Next we consider the lensing effect on the CMB polarizations, Q{n)ziziU (h). We are especially 
concerned with the rotationally invariant combinations, i.e., E- and i?-mode polarizations 
[30]; 

[E ± iB]L,M = [E± iB]L,M + Jdn ±23^2,m(^) [V0(n) + {-kV)w{n)] V(Q ± iU){n) 
= [E± ^B],,M + E E ± ^B],^^,{-l)'' (_^^ ^ O 



±2S^^ ll^ lliXi^rn , (2.9) 



X 



with the quantities, ±2^^^^, and ±25^^^/, defined by 



(-1)'' (_M L i) ^'Ke/ = I dn±2ylMi^)[Voyi,mm ■ [V±23^,'v(n)] , (2.10) 
(-1)*' (_M L m') ^'^le/ = I dn±2ylMi^)[{*V)oyiA'^)] ■ [V±23^,'v(^)]- (2.11) 
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The quantity -|-23^^' m' (^) denotes the spin-ib2 spherical harmonics. Similar to the case of 
temperature anisotropies, the quantities, ±2^^^^/ and ±25^^^/, are written as [30] 



(2^ + l)(2/ + l)(2L + l) 
167r 



+ l)+/(/ + l) -L(L + 1)] 



±2 =f2 , 



±2^ 



LIX' 



+ 1 T2 / L 



T2 / L 



' / + 1±2 V±2 -1 1t2; \//±2V±21-1t2 
Eq.(2.9) is rewritten in the separable form for E- and i3-mode polarizations: 

^-^ \ —M m m 

Bl,M = Bl M + } ' ( — 1)^ I ' 

where we define 



cx I cx 



qx 



2i 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



Note again that, for an even integer of L + ^ + /, the coefficients ©5™ and qS'^ vanish. On 
the other hand, the quantities ^iS*^ and e'^^ vanish when L+l+l' is an odd integer. Similar 
to the case of temperature, these properties come from the fact that i?-mode polarization 
and scalar lensing potential are even parity, while i?-mode polarization and pseudo-scalar 
lensing potential are odd parity. 



3 Reconstruction of deflection angle in the presence of curl mode 

In this section, we present a reconstruction method for (pi^m and von^rn^ based on the quadratic 
statistics (e.g., Refs. [21, 23-25]). We frequently use the formulas for Wigner-3j symbols 
summarized in appendix A. In what follows, the lensed temperature or polarizations, i.e., 0, 
E or are symbolically denoted by X (and Y\ 

3.1 Full-sky formalism 

3.1.1 Lensing field as quadratic statistics 
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Table 2. The functional forms of //^^ ^, and f^L l' appeared in Eq.(3.1). The label "even" and 
"odd" indicate that the function are non-zero only when ^ + i + i' is even or odd, respectively. 



XY 


h,L,L' 


h,L,L' 


ee 




o^ZiM^f^ - ^^^L^i^LCf" (odd). 




oKiM^L'"" + ^^t',,LCf'' (even) 


o^ZiM^L'"" - ^^l'al^l'' (odd) 


QB 


- e^J, ,,^Cf ^ (odd) 


Q^I'AL^l'' (even) 


EE 


^4,,L'CF + e5,%,,Cf ^ (even) 




EB 


-e<,L'^r-e5,^^,,Cr (odd) 


-e5- ,,Cr+e5^,,,^Cf^ (even) 


BB 







To get insight into the reconstruction of the lensing potentials, we consider an idealistic 
situation; we can take the ensemble average over primary CMB anisotropics alone, under a 
given realization of the lensing potentials. Hereafter, we denote this average by (• • •)cmB' 
distinguish it from the usual meaning of the ensemble average, (•••). Under the situation, 
in the correlation of lensed CMB anisotropics, {Xl^m^l' Af')cMB) the lensing potentials are 
included in the non-zero off-diagonal terms (L 7^ L', M 7^ —M'). This is because the 
lensing effect causes a non-trivial mode-coupling between the primary CMB anisotropics and 
lensing potentials (see Eqs.(2.3), (2.14) and (2.15)), and the lensed CMB anisotropics are not 
statistically isotropic for a given realization of the lensing potentials. Thus, it is possible to 
reconstruct the lensing potentials by extracting the off-diagonal terms of {Xl^mYl' m')cmb- 

To see this, let us write the correlation of the lensed CMB anisotropics in the harmonic 
space. With Eqs.(2.3), (2.14) and (2.15), we obtain 

with the coefficients, ff'!:^.'^^ and fT'}''^'i^\ summarized in Table 2. The conditions, "even" 
and "odd" , in Table 2 come from the parity symmetry in the lensing potentials and primary 
CMB anisotropics. To extract the off-diagonal terms and find the solutions for (pi^m and 
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wirn: we multiply 



in both sides of Eq.(3.1). Note that the multipoles, L and L' , are chosen so that fp^^^/^ 7^ 0. 

Then, summing up the equation over M and M' , and using the formulas, Eqs.(A.3) and 
(A.4), we find ^ 



£,L,L' M,M' 



CMB 



Notice that the term involving w in Eq.(3.1) vanishes. This is because, for all i, L and L' , 
the parity symmetry of ff'j^^J^^ and f^^^/^^ (Table 2) leads to 

/^'f/'/Slf'-O. (3.4) 

Similarly, following the procedure described in Eq.(3.2) below, but replacing f^^^"^ with 
f^^'^^\ the solution for tu^^m is obtained, and the result is 

■'l,L,L' M,M' 

The above equations, (3.3) and (3.5), can not be used for a definition of the estimator of 
lensing potentials, because these equations include the ensemble average over the primary 
CMB anisotropies alone, (• • •)cmb- But, the above equations imply that, by summing the 
quadratic combination of lensed fields over multipoles appropriately, it is possible to sep- 
arately construct the estimators for the scalar and pseudo-scalar lensing potentials, 4> and 
zu. 

3.1.2 Estimator 

Based on Eqs.(3.3) and (3.5), we first naively define the estimator for the lensing potential 
X (= (/) or w) as follows: 

■'e,L,L' M,M' 

where the subscript, a, means a pair of two CMB maps, e.g., a = G0 or EB. The multipoles, 



L and L', are chosen so that 7^ 0. With Eqs.(3.3) and (3.5), the estimator is rewritten 



as 



Eq.(3.6) = x,,„ + n^'^"^^,, (3.7) 



^ In deriving Eq.(3.3), we have ignored the zero-mode (Cq^), arising from the first term in Eq.(3.1) 
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where the quantity, n^,'^^\ r'j is given by 

■'£,L,L' A/, A/' 

Note that the above equation can be expressed without the quantity, (• • •)ciy[B- instance, 
in the case using the temperature anisotropics alone (i.e., a = 00), the above equation is 
rewritten as 



M 



e,L,L' MM' 



{^-^^^'' [-M' L' M'') ''^'l'al"^l",m"Ql,m 



Kml'm' L",M" 



+ (-1)*' (_M i' M'') o<,,^"0L",M"e^',M'}4',^'} ■ (3.9) 

The above estimator (3.6) suffers from several drawbacks in practical application to 
observation. At first, we should know about the primary CMB angular power spectra included 
in f^'^^i a priori. Another problem is that the estimator includes the contribution from the 

term n^'^^j^ which leads to a noisy reconstruction of the lensing potentials. Nevertheless, for 
the former point, the primary CMB angular power spectrum can be theoretically inferred if we 
know a set of fiducial cosmological parameters from other observations. On the other hand, 
for the latter point, we redefine the estimator for ^ and vj by introducing a weight function, 
-^f '/"rl, in order to reduce the contribution from n'^'''"^ ^i. Summing up all combination of L 

t,L,L t,m,L,L ^ 

and L , we write the estimator of the lensing potentials as 



L,L' MM 



The functional form of the weight function is determined so that the noise contribution is 
minimized. 

In what follows, we determine the functional form of the weight function so that the 
estimator is unbiased and the noise term is minimized. Eq.(3.10) can be recast as 

J 



= E (3.11) 

where the inner product [a^, ]^.^^ for arbitrary two quantities, a^'i^l, and 6^ '^^'1, is defined 

by 
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The quantity, n^'^\ is defined as 

L,L' 

Eq.(3.11) implies that the estimator would be an unbiased estimator if we impose the follow- 
ing condition: 

[^^r']?^ = '^.y• (3.14) 

Mathematically, this is equivalent to {x^f^)cMB = xe^m- Also, we wish to suppress the noise 
contributions, n^'^\ imposing the following condition: 

<t^\')=0. (3.15) 



Let us determine the functional form of the weight function under the conditions, (3-14) 
and (3.15), with the Lagrange- multiplier method. The variance of n'f^^ is given by 

/in^") |2\ - ^ VfF^''^"^ )* 
/ ~ 2^ + 1 

L,L' 

X (^L^'Cf ""'cT + F^^;iTi{-iY^'^'' ej' e^,^) , (3.16) 

where the quantity, C^^, is the lensed angular power spectrum including the contributions 
from instrumental noise. The detailed calculation for the noise variance, (I'^^'^P) , is presented 
in appendix B. Then, Eq.(3.15) under the constraint (3.14) is equivalent to 

L,L' 

+ E A^([F^r']?^-<5..')} = o. (3.17) 



"^LL,L' ^ L,L' 



The quantities, and A^, are the Lagrange multiplier whose functional form is specified 
below. Eq.(3.17) leads to 

In the above, interchanging L and L', we also obtain 

{F^Pjc^^'^cr + (F;f^) )*(-i)^+^+^'cT^c7/ + y: K'flt^l = • (3-19) 



X =0,ZU 



Multiplying the factors C^^C^^ and -(-1)^+^+^ Cf ^C^^ with Eq.(3.18) and Eq.(3.19), 
respectively, the sum of Eqs.(3.18) and (3.19) gives 
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where we define 

/ f^Xa) \*rXXrYY ( -\\i+L+L' pXYnXY ( f^Ao') \* 
x,ia) _ ^h,L,L'> ^L' " [-^) ^ l' ^h,L' ,l' , . 

^e,L,L' nXXnYY rXXnYY _ (pXY rXY\2 ' ^ ' 

Substituting Eq.(3.20) into Eq.(3.14), we obtain 

-Y.^K'W"^r'P =^xx'- (3.22) 

x" 

From Eq.(3.4), we find 

[f\rV = ^x"x'\f\rV- (3.23) 

Combining the above equation with Eq.(3.22), we obtain the exphcit form of the Lagrange 
multiplier 

(A^)* = . (3.24) 

Then, from Eq.(3.20), we finally obtain the expression for the weight function: 

x,{a) 

X (a) ^ JwJ (3 25) 

Note that, with the explicit expression (3.25), the noise variance, N^'^°'\ given in Eq.(3.16) 
becomes 



A7a;,(a) ^ / |„a;,(a)|2 



2^+1 [/^,5^l(5') Z^^^l^L^L') 



L,L' 



( ^x,{a) RXXnYY , T^x,{a) , ^\l+L+L' nXY pXY\ 



1 



i(") 



(3.26) 



where we use the relations given in Eqs.(3.18) and (3.24). Thus, the weight function can be 
recast as 

With the weight function given above, the estimators defined by Eq.(3.10) become opti- 
mal, i.e., the noise contribution is minimized. Eq.(3.25) or Eq.(3.27) is one of the main results 
in this paper. Note that, if the curl mode is absent, tz7 = 0, the resultant form of the weight 
function for (p exactly coincides with the one obtained in Ref. [25]. The difference appears 
when the angular power spectrum of the pseudo-scalar lensing potential, Cf^, included in 
the lensed angular power spectrum becomes non-vanishing. Note again that, in practical 
case, to use the estimator, the angular power spectrum of primary CMB anisotropics should 
be a priori known (i.e., ff'jflr and f^^^i are given). 
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3.1.3 Optimal combination 

As discussed in Ref.[25], the noise contribution can be further suppressed by combining multi- 
ple observables. Summing up the whole possible combination of temperature and polarization 
anisotropics, the optimal combination of the minimum variance estimators are given by 

-^M^^.W^W = (3.28) 



The weight functions, W^'^°'\ are determined so that the estimator satisfies the unbiased 
condition {{x^iI^)cmb = X£^rn), and the variance of the noise contribution is minimum. The 
optimal combination of the minimum variance estimator is then determined by the same 
analogy as in Ref. [25], and the result is 

.(c) _ Ar^,{c)V^r/ ;v-?^-U"'/5 



X 



a,(3 



where the variance, N^'^^\ is defined by 



The component of the matrix, {A?"|^}q,/3, is the covariance of n^'^^ and n^'j^'^ which is given 

by 



2£ + l A^.^' 



L.L' 



The derivation of the above equation is given in appendix B. 
3.2 Flat-sky limit 



The quadratic estimator for the curl mode has been empirically derived in previous 
work [32], based on the flat-sky approximation. Here, we show that our full-sky estimator 
can reproduce the flat-sky estimator of Ref.[32] (Eqs.(10)-(12) of Ref.[32]), in the flat-sky 
limit, £,L,L' < 1. 

Let us first rewrite the full-sky estimator (3.10) in Fourier space. In the flat-sky limit, 
we usually adopt the plane wave as a harmonic basis, and spin-0 quantity, A{n), such as, 0, 
and ro, is described in Fourier space as [34] 

Ai = ^ <fhe~'^-^A{h) , (3.32) 

where the two-dimensional vector, £, is given by {i cos (p£,i sin ipi). Similarly, in the flat-sky 
limit, E, and B modes are related to the Stokes parameters Q{'n.) and U{h) as [34] 

Ee±iBe = - [ (ffie^'^'^^-^'h-'^-^[Q ± iU]{h) . (3.33) 
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Table 3. Functional forms of /^'j, j^, and f ^ ^ £,/ in the flat-sky case. 



XY 


-fU^y) 

h,L,L' 


-rw,(XY) 

h,L,L' 


ee 


Cf Q£ • L + C®®£ • L' 


Cf^{-k£)-L + CfP{-k£)-L' 


e^ 


Cf^i ■ L cos 2ip^^^, + C^,^£ ■ L' 


Cf^{*£) ■ L cos 2ipL j^, + C®^(*£) • L' 


QB 




Cf^{-k£) ■ Lsm2tf^j^> 


EE 


[£ ■ LCl^ + I ■ L'Cf^p] cos 2(^^^^/ 


• LCf ^ + {*£) ■ L'Cf^p] cos 2(^^^^/ 


EB 


[I ■ LC^^ - e ■ L'C^P] sm2ip^^^, 


• LCf ^ - • Z'Cf,^] sin2(^^^^/ 


BB 


[£ ■ LC^^ + £ ■ L'ClP\ cos 2(^i^/ 


• LCf ^ + • L'C^,B] cos 2^5^^^/ 



Eqs.(3.32) and (3.33) imply that the Fourier coefficients (= yl^, ii^^ and Bi) are related 
to the harmonic coefficients, Zi^m (= ■^t,m-, Ei^m and Bg^m), through [34]: 



47r 



2£ + 



-^i-'-Z,,^e^-'^^ (3.34) 



2^ + 1.^ /-(/(/j^ 



47r / 27r 



^g-*m^,^^_ (3.35) 



Using the above equations (3.34) and (3.35), and combining Eq.(3.27), the full-sky estimator 
(3.10) in Fourier space is re-expressed as 



- / ^ / -^^t,L,L'^i 9,Il'^l^l' , (3.36) 

L,L' 



where we define 

_ / (2L + l)(2L^ + l) y/^ ^ (^^^(^ ^ ^ 

X g*"«/'<'e-«A'-''VLg-*A'/VL'^-m+M+M' ^ (3.37) 

To go further, we approximate the quantities, Ti^l,l'9£'l^i,' N^'^°'\ taking the flat- 
sky limit. To do this, we use the following relation valid under the flat-sky approximation. 



- 13 - 



^> 1 [34]: 



g±si(^,-^)gi^.n ^ ^^.^s (s = 0, 2) . (3.38) 



We also note that the delta function is given by [34]: 



e^^-^ . (3.39) 



(2vr)^ 



Using Eqs.(3.39), (3.38), (2.4), (2.5), (2.10), and (2.11), we find that, under the flat-sky 
approximation, 



U,L,L'9irr' - ^L+L' 



^L' yj£,L,L') ^L' l,L' ,L> 



cxxcYYc^xcyy _ ^c^ycxy^ 

/V^.(") ^ / /" / ^2W^ , 7^:(") -a=,(a) 

^Nf"\ (3.41) 

where the function, /^'^"^', is given in Table. 3 for each x and a. Note that the quantity, 

N^' , is the flat-sky counterpart of the minimum variance, . The detailed derivation 

of Eqs.(3.40) and (3.41) is given in appendix C. Then, Eq.(3.36) is rewritten as 

where we define the function, as 

^£,L,L'-^e 9^LL'- (3.43) 

Eq.(3.42) is the flat-sky counterpart of Eq.(3.10), which exactly coincides with that empiri- 
cally defined in Ref . [32] . Ref. [32] mentioned that their estimator does not satisfy the unbiased 
condition, and may detect some non-zero signals of vj even in the absence of the curl mode. 
But our result show that their estimator satisfies the condition, (x^"^)cmb = ^i, and becomes 
zero when the pseudo-scalar lensing potential vanishes ^ . 

Before closing this section, we also give the expression for the optimal combination of 
the flat-sky estimator used in the next section. With Eq.(3.35), the optimal combination 
(3.29) is rewritten as 



4^)=iVf(^)5;{(iV/)-lr'^4"^ (3-44) 

Q,/3 



^In fact, the integrand in Eq.(14) of Ref.[32] is an odd function in terms of the angle of £i, and the 
right-liand side of Eq.(14) vanishes. 
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Table 4. Experimental specifications for the PLANCK and ACTPol used in this paper. The 
quantity 61, is the beam size, and a^, represents the sensitivity of each channel to the temperature 
a^^T or polarizations ti^.p, depending on the power spectrum of temperature {X = O) or polarizations 
{X ~ E or B). The quantity v means a channel frequency. 



Experiment 


/sky 


V [GHz] 


Oij [arcmin] 


ay^T [/iK/pixel] 


au^P [^K /pixel] 


PLANCK [1] 


0.65 


30 


33 


4.4 


6.2 






44 


23 


6.5 


9.2 






70 


14 


9.8 


13.9 






100 


9.5 


6.8 


10.9 






143 


7.1 


6.0 


11.4 






217 


5.0 


13.1 


26.7 






353 


5.0 


40.1 


81.2 


ACTPol [3] 


0.1 


148 


1.4 


3.6 


5.0 



and the minimum variance is obtained from Eq.(3.30). Then, if we denote the flat-sky 

counterpart of N^''^'^'^ and 
flat-sky limit is described by 



counterpart oi and as N^' and , respectively, the optimal combination in the 



if = iV?'(=) 5]{(iV?)-n°'''4"^ , (3.45) 

a,/3 



with the variance, N^'^'^\ given by 



E{(^?)"'W- (3.46) 



^^e 13, p' 

The component of the matrix, {A7"^}„^^, is obtained by computing the flat-sky counterpart 
of Eq.(3.31) and the result is 



jx,(a,[S) _ f cfL f ^2t'x 

(27r)^ 



x(^S'crcr+^:i?:.crc^'n. (3.47) 



The detailed calculation of Eq.(3.47) is given in appendix C. 



4 Noise spectrum 

In this section, as a first step to estimate the feasibility to detect the curl mode based on the 
quadratic estimator, we compute the noise spectrum of the full-sky estimator, in the following 
cases; ACTPol combined with PLANCK (ACTPol-hPLANCK), and cosmic- variance limit 
(CV-limit). We also numerically evaluate the difference between the full- and flat-sky noise 
spectra. 
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Figure 1. The noise spectra of the scalar (left) and pseudo-scalar (right) lensing potentials for 
the lensing reconstruction from the temperature map alone, a = QQ (T; red lines), and from the 
temperature and polarization maps (T-l-P; green lines). We assume two cases; ACTPol combined 
with PLANCK (ACTPol+PLANCK; solid lines), and cosmic- variance hmit (CV-limit; dashed lines). 
We compute the noise spectra according to Eqs.(3.26) for T and Eq.(3.30) for T+P, with fmax = 3000, 
and take into account the effect of finite sky coverage. 
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Figure 2. The fractional difference of the noise spectrum between full- and flat-sky estimators (4.3), 
in the case with ACTPol+PLANCK. The left and right panels show the fractional difference for scalar 
and pseudo-scalar lensing potentials, respectively. 
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In the full-sky case, the noise spectrum is given by Eqs.(3.26) and (3.30). On the other 
hand, the noise in the flat-sky limit is obtained from Eqs.(3.41) and (3.46). To compute the 
noise spectrum, we assume that the lensing effect on the CMB comes only from the large- 
scale structure, and no source to produce the pseudo-scalar lensing potential is present. We 
use the modified version of FuturCMB code ^ [35]. In computing noise spectra, we further 
need the lensed, unlensed and instrumental noise angular power spectra. The lensed and 
unlensed power spectra are computed by CAMB [36] with the fiducial value of cosmological 
parameters described in section 1. Owing to the assumption, the angular power spectrum 
for the pseudo-scalar lensing potential, C^^, is set to zero. On the other hand, the angular 
power spectrum of the scalar lensing potential is obtained by integrating the matter power 
spectrum along the line-of-site, for which we adopt the fitting formula of the non-linear matter 
power spectrum given in Ref.[37]. The instrumental noise power spectra are given by [38] 



XX 



XX\-1 
Lv ) 



KfXX ^ ( 
■'^l,u — 7^ 6Xp 
V-'CMB/ 



2 



8 In 2 



(4.1) 



with TcMB = 2.7K being mean temperature of CMB. Here, the quantity 9y is the beam size, 
and represents the sensitivity of each channel to the temperature ay,T or polarizations 
Gy^p. The specific values for PLANCK and ACTPol are summarized in Table 4. Note that, 
for ACTPol-l-PLANCK, we assume that the survey region of ACTPol is entirely overlapped 
with that of PLANCK, and plot the following noise spectrum 



x,{a) 



ACTPol / ^PLANCK _ j^ACTPoh 

^sky v^sky •^sky / 



+ 



V^^£:ACTPol/ '^^^^iPLANCKV 



-1/2 

(4.2) 



where 1^"^^°^ = 0.1 and f^^^^^^ = 0.65 are the fractional sky-coverage of ACTPol and 
PLANCK, respectively, and the label, a, represents a pair of temperature maps, "00" or the 
optimal combination, "c" . In computing the noise of reconstruction within the survey region 
of ACTPol, ^I'xc^VoX'' lensed angular power spectra from PLANCK instead of 

ACTPol, aX, ^ < 700, in order to remedy a large uncertainty at large angular scales arising 
from the atmospheric temperature fluctuations. On the other hand, the noise, -^^p^anck' 
calculated with PLANCK experimental specification. For cosmic-variance limit, the recon- 
struction noise is computed with the instrumental noise being Mf''^ = 0. 

In Fig.l, we plot the expected noise spectrum in the full-sky case. The left and right 
panels show the noise spectra, Nf'^^^ and N^'^"'\ respectively. The resultant noise spectra for 
the pseudo-scalar lensing potential have amplitude comparable to those for the scalar lensing 
potential. In the cosmic-variance limit, the reconstruction noise is improved by more than 
an order of magnitude compared to the case with ACTPol-|-PLANCK. Notice that the noise 
spectra for the pseudo-scalar lensing potential is sensitive to the inclusion of polarizations 
compared to the estimator of the scalar lensing potential. In our calculation, the angular 
power spectrum of the pseudo-scalar lensing potential is set to zero. Although the primordial 
gravitational-waves or cosmic strings induces the pseudo-scalar lensing potential, as shown in 
the next section, the amplitude of is at least two orders of magnitude smaller than that 
of Cf^, as long as we consider the model parameters consistent with observations. Thus, the 
inclusion of the pseudo-scalar lensing potential would hardly change the result in Fig.l. 



http://lpsc.in2p3.fr/perotto/ 
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Figure 3. The angular power spectrum of the pseudo-scalar lensing potential from primordial 
gravitational-waves with the tensor-to-scalar ratio, r = 0.1 (top left), and cosmic strings with Gfi = 
10~®, P = 0.001 (top right). The error boxes in each figure show the expected variance of angular 
power spectrum from ACTPol combined with PLANCK (red) and cosmic- variance limit (green). The 
bottom two panels show the signal-to-noise ratio as a function of maximum multipole, for primordial 
gravitational-waves (bottom left) and cosmic strings with G/i = 10~*, P = 0.001 (green). 



In Fig. 2, to show the difference of the noise spectra between the flat- and full-sky 
estimators, we consider ACTPol-|-PLANCK and plot the following quantity as a function of 
£: 

Af''^ = . (4.3) 

The fractional difference becomes significant at large scales, and the difference becomes ^ 10% 
at 10 for both (j) and zu. The results do not sensitively depend on whether we include 
the polarization data for the reconstruction or not. Note that the fractional difference for (p 
is roughly consistent with Ref.[25]. These results show that the flat-sky estimator becomes 
invalid on large scales {l^ 10) and the full-sky lensing reconstruction is essential to extract 
information on primordial gravitational-waves and cosmic strings as discussed in the next 
section. 



5 Implications for primordial gravitational-waves and cosmic strings 

In this section, we illustrate the usefulness of the pseudo-scalar lensing potential as a diagnosis 
of the vector/tensor perturbations. Here we specifically focus on the pseudo-scalar lensing 
potential induced by two cases; primordial gravitational-waves produced during inflation, 
and cosmic strings (e.g., Refs.[39, 40]). 

The gravitational waves can produce the metric perturbations which have odd parity 
symmetry. This means that the lensing effect induced by the gravitational waves cause the 
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curl mode of deflection angle. There are several studies on the curl mode of deflections (or 
the pseudo-scalar lensing potential) induced by the primordial gravitational-waves ([30, 32, 
41, 42]), and the angular power spectrum of the pseudo-scalar lensing potential is given by 
[42] 



with 

AF^(fc,r?)= rfcd/3j,(fc,)^,,(-) 



(5.2) 

=k{rio-ri) 



kf] J 

The quantity, r/^,, denotes the conformal time at the last scattering surface, and the quantity 
r is the tensor-to-scalar ratio and we assume r = 0.1 which is close to the current upper 
bound [43]. 

On the other hand, extended objects such as cosmic strings produce not only the scalar 
perturbations but also the vector /tensor perturbations, and induce the pseudo-scalar lensing 
potential. For simplicity, we focus on the pseudo-scalar lensing potential only from vector 
perturbations. To compute the pseudo-scalar lensing potential, we consider a string network 
described by the velocity-dependent one-scale model [39, 40, 44-47], characterized by the 
intercommuting probability, P. The loop formation process is required for the energy loss of 
cosmic strings and the scaling solution. Also, this process is related to the number density 
of strings. Note that the constraint on P has an important role to distinguish between the 
cosmic strings as the topological defect {P = 1) and cosmic super-strings generated from 
the stringy inflation {P <^ 1) [48, 49]. We assume that the string segments are distributed 
randomly between the last scattering surface and observer, consistently with the scaling 
model, or, in other word, we neglect the correlation between two different segments [40]. 

If the string tension, G/i, and intercommuting probability, P, are given, the pseudo- 
scalar lensing potential from cosmic strings are computed as follows [29]: 



with [29] 




C^ = J(^^^^^^G^,)\^--^^ / -[A,-(fe,r/)]^ (5.3) 



( K 



5 
erf 



1/2 

k{x) 



=k{riQ-'q) 



(5.4) 



The quantity, H, is the Hubble expansion rate, and the quantities, ^ and v, are the correlation 
length and root-mean-square velocity, respectively, and determined from [39] 

2v^{l + — ^0.23P F =1, 5.5 

^ _ 2V2 l-%v^ 



The detailed description of model for the pseudo-scalar lensing potential from cosmic strings 
is described in the forthcoming paper [29]. 

In the left two panels of Fig. 3, the angular power spectrum from primordial gravitational- 
waves and the signal-to-noise ratio are shown. The plotted errors are estimated from 

AQ= ^^ -^ , (5.7) 

V(£ + l/2)/,kyA£ 
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where A£ is the size of multipole bin, and we set A£ = (z + 1)^ — for i-th bin, just for 
ihustration. For ACTPol combined with PLANCK, we evaluate the errors as 

= {(ACf:ACTPol) + (AQ:PLANCK , (5.8) 

where the errors arising from ACTPol survey region, AC^^^rpp^j, are computed according to 
Eq.(5.7) with /sky = 0.1 and = A^^fACTPoi- Similarly, the errors from PLANCK survey 
area, AC^j^^p^^j^, are obtained from Eq.(5.7) with /sky = 0.55 and = -^^lanck- 
the other hand, for the cosmic-variance limit, we compute the errors with /sky = 1-0 and 
the instrumental noise power spectra being zero. The signal-to-noise ratio for angular power 
spectrum of the pseudo-scalar lensing potential is defined by 




(5.9) 



As is expected, it is hard to detect the signature of primordial gravitational-waves from 
lensing reconstruction. For ACTPol+PLANCK, the signal-to-noise ratio is less than 0.1. 
Even with CV-limit, the signal-to-noise ratio is ~ 2. For the tensor-to-scalar ratio below the 
current upper limit [43], the signal-to-noise ratio would be further worsen. This is true as long 
as we adopt the quadratic estimator. These conclusions are consistent with the one obtained 
in Ref.[32], where they discussed the detectability in the flat-sky limit. Note that, recently, 
Ref . [50] discussed the detectability of the primordial gravitational- waves using the odd-parity 
bipolar spherical harmonics. Although they consider the temperature alone, their estimator 
(Eqs.(12) and (13) of Ref.[50]) and the variance (Eqs.(14) and (15) of Ref.[50]) coincide with 
our estimator (3.10) and the variance (3.26), respectively. However, the estimated amplitude 
of the noise spectrum for curl mode, assuming the PLANCK experiment alone, is several 
orders of magnitude smaller than our estimate based on ACTPol+PLANCK (red lines plotted 
in Fig.l). Then, they claimed that the PLANCK would detect the primordial gravitational- 
waves even from the curl mode if the signal is comparable to a currently detectable level, 
r ~ 0.24 [43]. This apparently contradicts Ref. [32] and our result. The discrepancy of the 
results may come from the experimental specifications of PLANCK adopted in Ref. [50] which 
may differ from ours, but the reason is unclear. 

On the other hand, in the right panels of Fig. 3, we show the case of the cosmic strings, 
with G/j, = 10~^ and P = 0.001 [51]. Note that these values are consistent with the con- 
straint from the temperature angular power spectrum using the Gott-Kaiser-Stebbins (GKS) 
effect [52, 53] induced by the gravitational potential of a moving string [39]. Although the 
result depends on the parameters of comic strings, the signal-to-noise ratio is ~ 3 even for 
ACTPol+PLANCK case. In the cosmic variance limit, the signal-to-noise ratio becomes 
~ 30. Note that, we have ignored the tensor metric perturbations from cosmic strings which 
also induce the pseudo-scalar lensing potential, and the inclusion of the contributions from 
the tensor perturbations would further increase the signal-to-noise ratio. The GKS effect 
observed via temperature map would have larger signal-to-noise ratio than the pseudo-scalar 
lensing potential [39]. However, the temperature anisotropics at small angular scales are 
usually dominated by the contributions from point sources and the Sunyaev-Zel'dovich (SZ) 
effect [54]. In this respect, the reconstruction of pseudo-scalar lensing potential is useful 
to check systematics and biases in the derived constraints on cosmic strings from GKS ef- 
fect. Note finally that the information on large scales would be important for detecting the 
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pseudo-scalar lensing potential from cosmic strings, and the full-sky formalism for lensing 
reconstruction is indispensable. 

6 Summary and discussion 

In this paper, we presented a full-sky algorithm for reconstructing the lensing potential of 
scalar (gradient mode) and pseudo-scalar (curl mode) components. We defined the estimator 
as a quadratic combination of observed anisotropies, and introduced the weight function to 
reduce the noise contribution (see Eq.(3.10)). The resultant form of the weight function which 
minimizes the noise contribution is given by Eq.(3.27) with Eqs.(3.26), and (3.21). Thanks 
to the different parity symmetry between scalar and pseudo-scalar lensing potentials, the 
gradient and curl modes can be separately reconstructed. Note that the quantities used to 
reconstruct the lensing potentials are summarized in Table. 1. In the flat-sky limit, we showed 
that the estimator reduces to the one empirically defined in Ref.[32]. We explicitly evalu- 
ated the noise spectra, and showed that the noise contribution for the pseudo-scalar lensing 
potential is comparable to that for the gradient mode. Further, prospects for reconstructing 
the curl mode is discussed, and signal-to-noise ratio for the pseudo-scalar lensing potential 
is computed, especially focusing on primordial gravitational-waves and cosmic strings. 

In this paper, we specifically focused on the lensing reconstruction based on the quadratic 
estimator proposed in Ref.[25]. On the other hand, in Ref.[26], for experiments sensitive to 
-B-mode polarization, it would be possible to improve the precision of the lensing poten- 
tial with the estimator based on the likelihood analysis. This may be also true for the 
pseudo-scalar lensing potential. Extending the estimator of Ref.[26] to the full-sky case, 
the signal-to-noise ratio estimated in section 5 would be improved, and upcoming or next 
generation experiments may detect the pseudo-scalar lensing potential even from primordial 
gravitational- waves. This will be investigated in our future work. 

Throughout this paper, we have assumed several idealizations, i.e., the higher-oder 
terms of the deflection angle are negligible, and observed CMB maps are given on the full 
sky without foregrounds and the inhomogeneous noise. However, at small scales, the lowest- 
order approximation (e.g., ignoring the higher-order terms in Eq.(2.2) for the temperature 
anisotropies) would not be valid for high resolution experiments [55]. The higher-order terms 
of deflection angle produce additional contributions in the lensed anisotropies, and the esti- 
mated scalar and pseudo-scalar lensing potential include the contributions from the higher- 
order terms of deflection angle [56]. Also, the masking effect [27, 57, 58], foreground contam- 
inations from point sources and thermal/kinematic SZ effect [10, 55], and the inhomogeneous 
noise [59] induce the additional non-zero off-diagonal terms in Eq.(3.1) and the estimated 
lensing potentials would be biased. Thus, in practical cases, to reduce the systematic bias in 
the estimated lensing potentials, an accurate treatment of these practical problems for the 
lensing reconstruction is required and worth investigating. 
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A Useful formula 

In this appendix, following Refs.[60] and [61], we summarize the formulas used in this paper. 

A.l Symmetry of Wigner-3j symbols 

Symmetric properties of the Wigner-3j symbols are described by 

nil "T-s J \"^2 "^3 n^i J \Tn-i mi m2 J \ms m2 mi y ' 

1 h ^ f-^Y^+£2+ia f ^1 ^2 ^3 A /^ 2) 

mi m2 T7i3 J y—mi —m2 —m^ J ' 

A. 2 Summation of Wigner-3j symbols 

Throughout the paper, we frequently use the following property of the Wigner-3j symbols: 



M 



^ I i L L'\ f i' L L'\ _ 1 ^ 

-m. M M' -m' M M' " WTl^^'^'V ■ ^^-^i 



m M M J \ -m M M J 21+1 

M,M' ^ 

A.3 Relations between Wigner-3j symbols and spherical harmonics 

Let us define 

^^rnV,Af ^ I dhsylMmVoye,m{n)] " [V .3^,.,^K^)] , (A.5) 

^^^^X,M ^ I dn,yi,,{h)[{*V)oyi,m{n)] . [Vsye'^m'in)]. (A.6) 
These two integrals are related to the Wigner-3j symbols as 



^m,m',M _ . /(2Ll + 1)(2L2 + 1)(2L3 + I) f ll h h \ fh h ^3 



^ V 167r \mi m2 in^ J \S —s^ 

X [-{ii + s)(4 - s + 1) + £2(^2 + 1) + (4 - s)(^3 + s + 1)] , (A.7) 



^m,m',M I {2il + I) {2£2 + I) {2h + I) f h I2 ^3 



L V 167r V"^i "^2 ms^ 



^J^2{^2 + lm+sm-s-l) 
^jm2 + mz-sm+s+i) 



s -1 -s + 1 



s 1 -s-1 



(A. 



B Noise covariance 



Here we derive Eqs.(3.16) and (3.31). Since Eq.(3.16) is obtained by setting /3 = a, in what 
follows, we derive the expression for the noise covariance given in Eq.(3.31). 
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B.l Relation between noise and estimator covariance 



First, we rewrite the noise variance using tiie variance of tiie estimator. Since n^'^^ 



i^ — Xi^rn, tlie noise variance is rewritten as 

(K:i"V<'r ) = (i-tlH'l) - - {{-tlr-^,m) + cr . (b.i) 

Note here that 

{{Xi,m)*^L,M^L',M') 

= EE E E 

x' l',m' L",M" Z=0,E,B 
X {(-1)''' (_^^/ ^ ^^//) 4;S'((x,,„)*Z^.,Jv/''4^„^'^^-,M) 

where we use Eqs.(2.3), (2.14) and (2.15), and introduce the coefficients, ■s^'^/^T^ defined by 

x,(ee) _ x,{EE) _ x,{BB) _ x,{EB) _ _ x,{BE) _ 

(B.3) 

and s^^f^i'^ = for the other combinations of XY . Assuming that the lensing potentials and 
primary CMB anisotropies are a random Gaussian field, and the correlations between the 
lensing potentials and primary CMB anisotropies are negligible, the above equation (B.2) 
reduces to 



'{x^^^rX^MYL'M) = i-^r M m) f'tZ'^Cr . (B.4) 

With Eq.(A.4), this leads to the following equation: 

>.™)'4S) = E ^l:^ E, (_i M m) (-1 M m) f^SP^r = cr ^ (b^s, 



L,L' M,M' 

As a result, Eq.(B.l) becomes 

^) = {{xtlfifl) - cr . (B.6) 

B.2 Estimator covariance 



Next we compute the estimator covariance, ( {x^"'^)* x^f^) ■ From Eq.(3.10), the covariance 



is given by 



.4a)y4i3)\^ \^ Sr^ \^ f ^ Li L[\ f i L2 

^^£,n,) ^e,m/ [-rn Ml M'J [-m M2 M' 

Li,L\L2,L'^Mi,Mi_' M2,M'2 

X (^S'/^ti^ (^2.,a/.^A,m;^^.,m.H^^^,mO . (B.7) 
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To compute the right-hand side of Eq.(B.7), we need to evaluate the four-point correlation 
of the observed CMB anisotropies: 



^^Lr,Mi^L'^,M\^L2,M2WL'^^M'^j ■ (B.8) 

In general, the four-point correlation (B.8) is described by 

(^Li ,Mi ,M[ ^L'^,M'2 ) 

The first term is defined by 

(^Li ,Mi yl'^ ,M[ ^L2,M2 ) Q 

~'~'Li '-Li' ^(l,l').{2,2') +'"^1 '"Li' ^(l,2),(l',2') + '"Li '^Li' ^(l,2'),(l',2) ' [D.W) 

where the quantity C is the observed angular power spectrum, and A^-^ ^2 2')' ^(1 2) (1' 2') 
and A(^i 2') (1' 2) denote 

^(l,l')>(2,2') = huLi'^L2,L'.,^hh-M\h.h-M'^ ' (B-H) 
^(l,2),(l',2') = ^LuL2h[,L'2^MuM2^M[,M2 ' (B.12) 
^(1,2'),(1',2) = hi,L'2h[,L2^Mi,M2^M[,M2 ' (B.13) 
The other terms included in the four-point correlation are represented by the second term 
in Eq.(B.9). If the quantities, Xl-^^^Mi,, Yl' m' > ^L2,M2 ^^'i ^l' m' > random Gaussian 
fields, the second term vanishes. From Eq.(B.9), the covariance is decomposed into two parts: 



where we define the Gaussian and connected parts, G^'^'^^ and C^'^'^\ as 



QX,[a,P) ^ ^ ^ ^ ^ i Li L['^ ^ i L2 L'2 

Li,L[ L2,L'^ Mi,M[ M2,M'2 



m Ml M\J \-m M2 M'^^ 

^x,{a,i3) ^ sr^ f i Li L[ \ / i L2 L'2 

W,m \-m Ml M'J \-m M2 M' 

Li,L[L2,L'2Mi,M'tM2,M'2 

Let us first compute the Gaussian part, G^'^'^\ Substituting Eq.(B.lO) into Eq.(B.15), 
the Gaussian part of the covariance is given by 



^ f ^ Li L\\ f I L2 L'2\ 

^i,ra ^ , ^ , V-m Ml M[J \-m M2 My 



Li,L'i L2,L'2 Mj,M'i M2,M'2 

X ^^i,L,,L[> ^IMAX^^ ^(I'l ) = (2.2') 

+ Cf,^C^^^A(i 2),(l',2') +<^Li^C^/^{l,2'),(l',2)} • (B-17) 
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Using Eq.(A.3), the term proportional to A^-^ ^^2') gives 61^, and we neglect this term to 
consider £ > 0. Then, Eq.(B.17) becomes 

QX,(a,l3) = ^ ^ ^ X] ^ ^ Li L[ \ f i L2 L2 
Li,L[ L2,L'2 Mi,M'i M2,M'2 



m Ml M\) \-m M2 M'^ 



px,{li) nXZ(=iYW ( ^ Li L[ \ / I Li L\ 
IMA ^1 \-m Ml M[J \-m Mi M[ 

.x,W pXWnYZ ( ^ Li L\ \ ( I L\ Li 



Using Eq.(A.l), the above equation reduces to 



E 



1 Li L[ \ f i Li L\ 
-m Ml M'J \-m Mi M[ 



1 ^z?^-W 



21 + 1 iM,L[^ 

( j^x,{l3) pXZpYW j_ T^x,il3) , -.x£+Li+L[pXWpYZ\ m 1 

Note that we use (A. 4) from the first to the last equation. 

Next we compute the connected part of the covariance (B.16). Even if the primordial 
CMB temperature and polarization anisotropies are random Gaussian fields, the connected 
part of the four-point correlations is arising from the secondary effects such as the weak 
lensing. Assuming that the connected part is arising from the lensing effect, the connected 
part of the four-point correlations is given as 



, Ml ^l; ,m; ^l2,m2 ^a M'2 1 c 

- 1^^^ Ml M'J \-m" M2 M'J mAJ V,L2,l^^^" ' ^^'^"^ 



Other terms included in the connected part of four-point correlation, such as the non- 
Gaussian terms introduced in Ref.[62], induce additional noise term in Eq.(3.31), but the 
terms may be an order of mag nitude smaller than G^'^"'''^ [56,62,63]. Substituting Eq.(B. 20) 
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into Eq.(B.16), the connected part of the covariance is rewritten as 

^x,(a,/3) _ \^ f ^ Li L[ \ f £ L2 L'^ 



^x,(a,(i} ^ sr^ sr^ sr^ \^ f ^1 \ 1 ^ -^2 ^2 

e,m \-m Ml M'J \-m M2 M' 

Li,L[L2,L'2Mi,M'iM2,M'2 

x*pX,(/3) v;^ / Li L[\ f i" L2 4 

£M,L'J tMA ,f^„ \-m" Ml M[J \-m" M2 M'^ 

I ,m 

X 

With Eq.(A.4), the above equation is rewritten as 

e,m - 1^ 1^ 1^ 1^ 2£ + l 2i + l 



Li,L[L2,L'^e",m" X 

y N*pX,(/3) (j:x',(a) „x',(/3) ^x' ,x' 

X 

Using Eq.(3.14), we obtain 

Finally, substituting the resultant form of the Gaussian (B.19) and connected parts (B.23) 
into Eq.(B.6), we obtain Eq.(3.31). 

C Flat-sky limit 

In this appendix, using Eqs.(3.34)-(3.38), we derive Eqs.(3.40), (3.41) and the noise cross- 
spectrum in the flat-sky limit (3.47). 

We first consider the expression for 7g £, l'9^'l^l' (3.40) in the flat-sky limit. From 
Eq.(3.21), we obtain 

' t,L,L'9i^L,L' ~ nXXpYY nXXpYY _ (nXYpXY\2 ' ^ ' 

Note that, we use 

W' = (-l)'+^+^'V.^. (C.2) 
Then, we need to compute the quantity, (7^ ^ l')* f^ffl'- shown in Table. 2, this quantity 

includes {T^^lM^ ^'^l' e L 0't,L,L'y ''■^1 1 l' "^"^^^'^ s = or ±2. From Eqs.(2.4), (2.5), 
(2.10), and (2.11), we obtain 



m,M,M 



X ^m-M-M- 



J (fnsyl, _^f{n)[Voy£,-m{n)] 0^ [V ,3^L,M(n)] , 

(C.3) 
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where, for arbitrary two vectors, a and b, we define the products, 0^ and 0^, as 
a b = a ■ b, a 0^ b = {-ka) ■ b = — {-kb) ■ a . 



(C.4) 



Note here that, the right-hand side of Eqs.(2.4) and (2.10) is a real number, since the left- 
had side of these equations is a real number. Similarly, the right-hand sides of Eqs.(2.5) and 
(2.11) is a purely imaginary number. Thus, Eq.(C.3) is rewritten as 



/(2L + 1)(2L' + 1) 



(r,,i,r)%5^',,,L - (^4^(2£ + i)(LL') 



1/2 



E 

m,MM' 



X I 



m-M-M' / 72 - 



I d'hsyL'^^M'inWsyiMn)] 0x [Vo3^;,_^(n) 



/ ^(2L + l)(2L^ + l) 
V 2(2£ + l)LL' 

m,AI,M' 



1/2 






X 



(27r)2 ^-"^'.A/' 

Using Eq.(3.38) and assuming I, L, L' ^ 1, the above equation reduces to 



(C.5) 



(C.6) 



From Eq.(3.39), we obtain 



1 (s = 0) 

cos2ip^j^i {s = ©) 
^- sin 299,^/ (•5 = 0) 



(C.7) 



where we denote ^i^i> = — ^l'- Similarly, the flat-sky counterpart of {Ti l i^')* s<S^ ^ j^' is 
obtained by interchanging L and L' in Eq.(C.7) if £ + L + L' is an even integer. If £ + L + L' 
is an odd integer, we further multiply it by minus sign. From Eq.(C.7), we can define the 
following quantity: 



-jx,(a) 

JllM°l+l' 



(C.8) 



The functional form of /^'^"^' is summarized in Table. 3 for each x and a. Substituting 
Eq.(C.8) into Eq.(C.l), we obtain the following expression: 



T- a;, (a) r- 



rxxpyy rxxpYY _ (nxv pXYvi 



L ^L' 



r —a;, (a) 



(C.9) 
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with the quantity, being 

_r,(a) ^ ^L' U e,L,L' ' ^ L' t,L' ,L> 

^l,L,L' qXXqYYqXXqYY _ (C^ycXY\2 ' V ■ j 

Next, we consider the flat-sky counterpart of Eq.(3.31), and show Eqs.(3.41) and (3.47). 
Using Eq.(A.4), the covariance (3.31) is rewritten as 



m,M,M' m'M",M"' 



21 + 1 \-mM M'l \-m' M" M'" 



^* j^xm ^XX'^YY' j^x,(l3) , -..e+L+L'^XY'^X'Y 

L,L' 

2^ 2^ 2£ + 1 [-m M M'l [-m M" M'" 

m,M,M'm',M",M"' ^ ^ ^ 

27r y 27r J 2tt 



X 
L,L' 

where, for arbitrary integers. Mi and M2, we use the following equation: 



SmuM2 = J ^e-*^^"^-^"^^'^. (C.12) 

With the quantity defined in Eq.(3.37), the expression for the covariance (C.ll) is rewritten 
as 

X (F-(J)* ^'C7' + FliP^i-lY^^^^'cr'c^J"^) . (C.13) 



Note here that, in the flat-sky limit, we can define the following quantity: 

h^L'^^l^^L - r,,L,L'Fli% . (C.14) 

This is because, from Eqs.(C.9), the right-hand side of the above equation is proportional 
to the delta function, 5£,_|„£,'„f. Using Eq.(C.14), and 60 = l/vr, and assuming £,L,L' 1, 
Eq.(C.13) becomes 

^XX'^YY' ^XY'^X'Y\ ,^ . 
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Note that, in the right-hand side, we can choose two-dimensional coordinate system for the 
variables of integration, L and L\ so that the integrand of ipi does not depend on ipi. Then, 
the above equation reduces to 

The noise in the flat-sky limit, N^'^'^\ is obtained by a = /3 in the above equation: 
Note that the quantity -F"^'^"]^' is described by 

^£,L,L'-^i ^^-^^^ 

Substituting Eq.(C.18) into Eq.(C.17), we obtain the expression for the noise spectrum in 
the flat-sky limit: 
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